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A METHOD FOR DETERMINING FIELD RADIATION LEVELS FOR A 

RADIATING DEVICE 
FIELD OF THE INVENTION 

The present invention relates to radiating 
5 devices such as antennas* 

BACKGROUND OF THE INVENTION 

For health and safety reasons radiation levels 
from antennas must be accurately monitored and controlled. 
Regulatory bodies require radio frequency wireless 
10 operators to ensure radiation levels comply with the 
relevant health and safety standards in place. This 
requires the collection of technical data relating to 
radiation levels and exposure limits around a relevant 
antenna. 

15 Existing techniques for radiation estimation 

approximate the characteristics of the antenna to a point 
source and use far field theory to calculate radiation 
levels for an antenna. 

The present invention is aimed at providing an 

20 improved method for determining radiation levels for 
radiating devices such as antennas. 
SUMMARY OF THE INVENTION 

According to one aspect of the present invention 
there is provided a method for determining field radiation 

25 levels for a radiating device comprising the steps of 
determining far field radiation characteristics of a 
radiating device, providing a model of the radiating 
device , which model approximates the determined far field 
radiation characteristics and determining a near field 

30 radiation characteristic from the model for at least one 
point in space. 

It is preferred that the model is chosen to 
approximate near field radiation characteristics. 

It is preferred that the method includes the step 

35 of determining a boundary between near field and far field * 
radiation. 

The method may include the step of determining 
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near field radiation density from the model. 

Preferably the method includes the step of 
determining near field radiation power density level over 
a plurality of positions in space. It is to be understood 
5 that space includes occupied (physical structure present) 
or unoccupied space and any particular area around the 
radiating device. 

It is preferred that the method includes the step 
of radiation pattern or gain characteristics of the 
10 radiating device from the two orthogonal far field 
radiation patterns. 

It is preferred that radiation pattern or gain 
characteristics are determined from documented data such 
as that available in handbooks. 
15 Preferably the method includes the step of 

determining 3dB beam width in the two orthogonal far field 
radiation patterns - 

It is preferred that the method includes 
determining physical characteristics of the device to 
20 determine the far field radiation characteristics. 

Preferably the method includes the step of 
providing a model including representing the device by a 
plurality of radiation sources. 

According to one embodiment of the invention the 
25 radiating device comprises a wire antenna. 

The method may include the step of providing a 
model including estimating the length and spacing of each 
wire element forming the wire antenna. 

Each radiation source may comprise one wire 

30 element. 

Preferably the method includes calculating mutual 
coupling between all wire elements of the radiating 
device . 

Preferably the method includes assembling an N by 
3 5 N impedance matrix for the radiation sources. 
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The method preferably includes calculating 
voltage for each wire element. 

The method may also include determining current 
in each wire element . 
5 Preferably the method includes multiplying an 

inverse impedance matrix by a column voltage vector (zero 
for parasitic elements and 1 volt for driven elements) . 

The method preferably includes assigning one or 
more Huygen's wavelet point sources to each wire element. 
10 The method may also include calculating magnitude 

and phase of each point source from current determined. 

The method preferably includes assuming a sin (6) 
dependence for each wire element ((9) being measured from 
the direction of the elements) and summing the 
15 contribution of each point source to each point in space 
within the near field radiation pattern. 

It is preferred that the method includes both 
near field and far field effects when calculating the 
contribution of each point source. 
20 The method preferably includes scaling field 

strengths determined at each point in space, by power 
supplied to the radiating device. 

The method preferably includes radiating devices 
which are Yagi-Uda, log periodic , single or phased arrays 
25 of monopoles, dipoles, rhombic antennas and other regular 
or irregular wire antennas . 

It is preferred that the method includes 
providing a single point source for each wire element with 
a length less than half the wavelength of radiation 
3 0 emitted from the radiating device. 

The number of point sources for a radiating device will 
preferably be a multiple of the wavelength divided by two, 
i.e. twice the length of the radiating device divided by 
the wavelength. 

35 Preferably the method includes providing scaling 

factors for miscellaneous effects such as weather, 
obstacles, other antennas, metallic structures and 
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dielectric structures or other factors which affect 
radiation characteristics. 

According to another embodiment of the present 
invention the radiating device is an aperture antenna. 
5 Preferably the method includes determining 

physical characteristics of the aperture antenna from the 
beam width characteristics. 

The physical characteristics include physical 
size, impedance, the size of the aperture and field 
10 distribution. 

The step of providing a model may include 
representing the aperture by at least one Huygen's wavelet 
source . 

Preferably the method includes the step of 
15 summing the contribution from each wavelet source at each 
point in space. 

It is preferred that the method includes the step 
of summing the contribution from each wavelet source over 
a three dimensional coordinate system in space, e.g. 

2 0 rectangular, circular or polar coordinate system. 

Preferably the contribution from each wavelet 
source determined includes power, voltage and current. 

It is preferred that power density level at each 
point in space is determined using known power density 
25 formulas such as those described in more detail 
hereinafter . 

According to another aspect of the present 
invention there is provided a method of estimating 
radiation power density of electromagnetic radiation 

3 0 comprising the steps of identifying a radiating device, 

representing the radiation device as a plurality of point 
sources which radiate electromagnetic radiation, 
estimating power density level at a plurality of positions 
in space for each point source and determining the total 
3 5 power density level at each position, by summing the 
contribution of each point source to the respective 
positions in space. 
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It is preferred that the method includes 
displaying the power density level for a plurality of 
positions . 

Preferably the method includes displaying the 
5 power density levels in graphical, tabulated, 
diagrammatical, pictorial or other form. 

The method may include choosing positions which 
fit into a two-dimensional or three dimensional coordinate 
system. 

10 Preferably point sources include any portion 

which requires at least two portions to represent the 
radiating device. 

Preferably the method includes summing the power 
density level determined at one position for all point 
15 sources representing the radiating device prior to 

estimating the power density level at another position. 

Alternatively the method includes estimating the 
power density level for the plurality of positions for one 
point source, storing the estimated power density levels 
20 then estimating and storing the power density levels for 
the plurality of positions for another point source and 
sximming stored power density levels for each point source 
at each position to calculate a resultant power density 
level at each position. 
25 It is preferred that the method includes 

identifying a plurality of radiating devices. 

The method may include representing each 
radiating device as a plurality of point sources. 

Preferably the method includes calculating the 
3 0 total power density level for each radiating device. 

The method may include calculating the total 
power density level for each point source at each position 
and summing power density levels calculated at each 
position. 

35 Preferably the distance between each point source 

is determined by the distance between points in space. 
Thus in a rectangular coordinate system the spacing 
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between grid points will determine the spacing between 
point sources, with closer distances between adjacent 
points in space resulting in more point sources for the 
radiating device. 
5 It is preferred that the electromagnetic 

radiation measured is radio frequency electromagnetic 
radiation. However other radiation is included. 

Preferably the method includes calculating far 
field and near field tapering characteristics for the 
10 radiating device. This may be accomplished using 
predetermined formula or field measurements. 

Preferably the error measured is radio frequency 

error. 

Preferably the method includes calculating the 
15 far field and near field tapering characteristics for each 
position. 

Preferably far field Pd formula is 

CJ+2A5 

Powerat Antenna * 10 10 



20 



Gd = Antenna gain with respect to dipole of 
analysis angle 

PaA = Power sent to antenna after lossy items 
Di = Distance from antenna. 
Preferably D far field distance = 

271 

2(th€O n ***** " n wfFprthwa nrrtu rt>\ 



25 if aperture > 1A 



16 times 



*i <=>i/i n^i^oations - define break point of 

if aperture > J A 

Pd paraNF « 41.3*Pd (para - parabolic aperture) 
Breakpoint = .16 Dff - parabolic aperture 
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= .25 Dff - rectangular aperture 

Pd Rec tNF = 

According to another embodiment of tlie present 
invention there is provided a method of determining field 
5 radiation levels for a radiating device comprising the 
steps of determining far field radiation characteristics 
of a radiating device, determining the boundary between 
near field and far field radiation, determining the 
displacement of a point in space relative to the closest 
10 point on the radiating device and calculating the power 
density level at the point in space. 

Preferably, the method of determining field 
radiation levels for a radiating device utilising the 
closest point method includes the step of determining 
15 characteristics of the radiating device in a similar 
manner to that used in connection with any one of the 
previous aspects or embodiments of the invention. 

It is preferred that the power density level is 
determined by the aforementioned formula for Pd. 

2 0 It is preferred that the power density formulas 

are modified according to modification factors affecting 
the gain of the radiating device and the degradation of 
radiation as a function of the displacement of the point 
in space from the radiating device - 
25 The method may include modeling the radiating 

device as a plurality of point sources. 

Preferably displacement is determined by 
determining X, Y, Z vectors in space. 

Preferably the displacement is determined using 

3 0 azimuth and elevation angles of the point in space. 

It is preferred that the method includes 
determining the orientation of the radiating device in 
space to determine down tilt. 

The method may include the step of determining 
3 5 any offset in the displacement as a result of down tilt of 
the radiating device. 

The method preferably includes the step of 
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determining if the point in space is outside the width 
plane or length plane or height of the radiating device. 

Preferably the method includes the step of 
calculating the effective reduction of antenna aperture as 
5 a result of the displacement of the point in space from 
the radiating device. 
BRIEF DESCRIPTION OF THE DRAWINGS 

A preferred embodiment of the present invention 
will now be described by way of example only with 
10 reference to the accompanying drawings: 

Figure 1 shows a schematic diagram of a method of 
determining power density of a point in space determined 
using a closest point algorithm; 

Figure 2 shows a schematic diagram of a method of 
15 performing power density of a point in space determined 
using a Hyugen's wavelet method for an aperture antenna; 

Figure 3 shows exposure limit boundaries of a 
test antenna using traditional modeling techniques; 

Figure 4 shows exposure limit boundaries for the 

2 0 same test antenna using a closest point algorithm 

technique; 

Figure 5 shows a power density plot for the test 
antenna exposure limit boundaries shown in Figure 3; and 

Figure 6 shows a power density plot for the test 
25 antenna modeled according to the closest point algorithm 
in accordance with the present invention. 

The best mode of performing the present invention 
will be described with reference to three different 
embodiments • 

3 0 Antennas typically fall into two main categories. 

One of these is wire antennas and the other is aperture 
antennas . 

A method for determining power density radiation 
levels for aperture antennas in accordance with the 
3 5 preferred embodiment of the* invention may incorporate a 
closest point algorithm technique or a multiple point 
source technique . 
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POWER DENSITY CALCUIATIONS 

As described in the Microwave Engineers Handbook, 
Volume 2, Artech House 1971 to Sadd, Theodore [1] , the 
simplest way of estimating the power density radiated from 
5 antennas is to apply the far field power density formula 
to a point source representation of an antenna. 

To achieve sufficient accuracy manufacturers far 
field gain patterns must be used. These exhibit the 
antenna far field gain characteristics for all directions 
10 (i.e. 0 to 360°) in the horizontal and vertical planes. 
The far field power density formula is given by 

PoweratAntenna * 10 10 
4*Di 2 



Where Pd is the estimated power density, Gd is the 
15 antennas gain with respect to a dipole at the analysis 
angle, PoweratAntenna is the power sent to the antenna 
after lossy items such as a signal feeders, Di is the 
distance from the antenna. Units for the formula are 
Watts per centimeter squared . 
2 0 The inventor has observed that this far field 

estimation grossly overestimates the power density in the 
near field area around the antenna. From observations and 
theory, according to one embodiment of the invention an 
algorithm has been developed for determining how far away 
25 the far field is from an antenna (and hence when the far 
field power density formula becomes accurate) . The 
algorithm is described in greater detail hereinafter. 

If the effective aperture of the antenna is less 
than 1 wavelength at the operating frequency then: 
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_ wavelength 

If the effective aperture is greater than or 
ecjual to 1 wavelength at the operating frecjuency, then: 



2*Aea 

D 



wavelength 

wherein Aea is the antenna's effective aperture. 

NEAR FIELD CALCULATIONS 

To compensate for the overestimating effect of 
just using far field calculations it is necessary to 

10 determine the reduction in gain in the near field due to 
the finite size of the aperture. From power density 
observations for uniform line sources and tapered 
illumination aperture antennas a general break point 
distance must be determined where calculations change 

15 modes between near field and far field. An associated 

taper method can be used to calculate near field radiation 
levels. The break point for parabolic antennas is 
determined as 0.16 times the far field distance. The 
taper method for the parabolic antenna is defined as: 

20 Pd porahtthcNF ^4\3*Pd 

where the Di used in the Pd calculations is just D far- 
field, irrespective of how close a distance to the 
antenna. The break point for rectangular aperture 
antennas is defined as 0.25 times the far field distance. 

25 The taper method for the rectangular antenna is defined 
as: 

where the Di 2 used in Pd is replaced with Di times the 
breakpoint distance. This effectively makes the power 
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density decay with a rate of 1/Di in the near field rather 
than 1/D± 2 (which occurs in the far field) . 

These taper methods can be seen to still 
overestimate in the near field, but to a much lesser 
5 extent. This small overestimation builds some tolerance 
into the calculations to cater for antennas that may have 
differing properties different to those exhibited. 
APERTURE ANTENNAS 

Before power density levels can be determined for 
10 an antenna it is necessary to identify key characteristics 
of the antenna. This enables the antenna to be accurately 
modeled. 

According to one embodiment the far field 
radiation pattern can be determined from the Fourier 

15 transform of the electric field across the aperture. When 
most of the radiation is concentrated in a highly 
directional beam, then close to the antenna there is no 
significant diffraction pattern when observed close to the 
main beam direction. That is radiation from any point 

20 source in the aperture travels almost the same distance as 
radiation from any other point source in the aperture when 
the observation point is confined to the non- shadow 
region- the aperture dimensions translated in the 
direction of the radiation, a simple inverse square law 

25 can be applied and a good approximation to the radiation 
levels is obtained. In the shadow, i.e. outside the 
aperture translated in the direction of the main beam, 
diffraction effects occur and the field is attenuated. 



3 0 falling on an irregular hole in an opaque material. When 
observed very close to the opaque material, the shape of , 
the light beam is the same as the irregular hole (the 
Fresnel region) . When observed far from the irregular 



The optical equivalent is a parallel ray of light 
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hole, the image is that of a circle (Fraunhofer region) . 

In order to calculate the near and far field 
radiation density given only two orthogonal radiation 
patterns and the power supplied to the aperture antenna, 
the mechanical structure must be determined from these 
radiation patterns. Once an approximate aperture 
distribution is known the antenna is modeled by point 
radiators, ("Huygen's wavelet sources) across the 
aperture . 

According to the closest point algorithm method a 
full three-dimensional analysis of the antenna is achieved 
over the antennas length and width. The closest point 
algorithm is set out as follows: 

1. Find the displacement of the analysis point 
and the centre of the radiating element. 

2. Calculate the effective reduction of 
antenna aperture as a result of the 
displacement . 

3. Find the displacement vectors x, y and z 

0 dimension values between the centre of the 

antenna and the point of analysis. 

4. Rotate the displacement vector by the 
antenna's mechanical downtilt in three- 
dimensional space (to calculate the 

5 analysis with respect to the antenna face) . 

5. Find the Azimuth angle (the angle in the XY 
plane from the centre of the antenna's 
aperture face to the analysis point in the 
X-Y plane) ; 

0 6. Find the Elevation angle (the angle in the 

2 plane from the centre o£ the antenna's 
aperture face to the analysis point in the 
Z plane) ; 
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analysis height going along the antenna 
height plane (obtain through the rotated 
displacement's z dimension); 

8. From the rotated displacement , obtain the 
analysis distance along the antenna width 
plane (obtained through the rotated 
displacement's xy plane vector's component 
along the antenna face 

9. If the offset (analysis displacement height 
+ electrical downtilt compensation is 
inside +/-1/2 the effective antenna height, 
then the new source Z position is set to 
the antenna's centre height + the Z 
component of the downtilted offset 
Else 

The new Z position is set to the sign of 
the offset *l/2 the effective antenna 
height; 

10. If the antenna width plane offset is inside 
=/-l/2 the effective antenna width, then 
the new source X and Y positions are set to 
the antenna's centre X and Y positions + 
the width offset along the antenna face + X 
and Y components of the downtilted 
elevation offset , 
Else, 

The new X and Y positions are set to the 
sign of the offset *l/2 the effective 
antenna width. 

0 The new antenna source positions are then used in 

the calculation at this analysis point. 

In the summary the closest point algorithm 
determines that for any point in space either within or 
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outside the near field/far field radiation boundary, the 
power density level at the point in space is measured 
relative to only one point on the aperture antenna. This 
point is the closest point according to the displacement 
5 data which is determined from measuring the X, Y # Z 

coordinates of the point in space as well as the elevation 
angle and azimuth angle. For each point in space the 
radiating point of the aperture is merely shifted to the 
most appropriate position to become the closest point on 
10 the antenna to the point in space. 

Referring to Figure 1 for points in space lying in 
the Far field 10 power density calculations can be made by 
simply using the previously highlight power density 
formula. However for points in the Near field 12 inside 
15 the boundary 11 the taper function must be applied to 
accurately determine power density level. This taper 
function linearly reduces the effective aperture size in 
length and width as the distance from the antenna 
increases, until the Far field distance is reached, at 
20 which point the closest point algorithm merges back into a 
point source analysis. The Near field calculations are 
based on the modifying factors incorporated into the power 
density Near field formulas hereinbefore described. 

The closest point algorithm also takes account of 
25 the orientation of the antenna including its tilt. 

For any point 13 in space which is outside the 
face of the antenna aperture, the closest point source 
would be on the perimeter of the aperture 14. 

As an alternative to the closest point algorithm 
3 0 method power density levels can be determined using an 
accumulated point source method. 

Thus as shown in Figure 2 a point in space 15 has 
a power density level which is determined by the 
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accumulated effect of point sources 16 on the antenna 
aperture 17 . 

As with, the closest point algorithm method the 
aperture antenna must first be modeled from the two 
5 orthogonal far field radiation patterns. The size of the 
antenna aperture is then calculated in conjunction with 
data on the physical size of the structure. 

Once the size and shape of the aperture has been 
determined the aperture is then represented by a number of 
10 Huygen's wavelet sources. The number of wavelet sources 
may be determined by using a single point source for each 
element with a length less than half a wavelength. For 
elements longer than this the number of point sources 
increases proportionately with the number of half 
15 wavelengths in the element. 

As an example the Yagi-Uda style antenna is 
discussed. The 3dB beam width in the two orthogonal 
radiation patterns are used to gain an estimate of the 
number of elements in the Yagi-Uda antenna. A standard 
20 configuration is then used to represent the antenna. 

The sum of the contributions from each point 
source to every point in space in the vicinity of the 
antenna is used to determined the Near and Far field 
intensities. Thus as shown in Figure 2 the contributions 
25 from each of the point sources 16 are combined to provide 
a power density level at point 15 . 

As with the closest point algorithm the field 
strengths are scaled in accordance with the power supplied 
to the antenna. 
3 0 EXPOSURE LIMITS COMPARISON 

The following diagrams display the exposure limit 
boundaries (using 2 and 10 W/m 2 for the red and yellow 
zones respectively) of the existing and new calculation 
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techniques for plan (horizontal) and elevation (vertical) 
boresignts of the test antenna, the Argue CTA610D-R. Note 
that in the traditional full modeling solution the antenna 
is identified as a point source located at its phase 
5 centre. See Figure 3 

POWER DENSITY COMPARISON 

The following diagrams represent the power density- 
graphical representations of the existing and new 
algorithms for the test antennas, of the type Decibel 

10 DB580y (omnidirectional) . Note that in the traditional 
full modeling solution, the antenna is identified as a 
point source located at its phase centre. See Figure 4 
WIRE ANTENNAS 

The Huygen's wavelet method can also be applied to 

15 wire antennas. Such antennas can be classified as Yagi- 
Uda, Log Periodic , Phased Arrays of Monopoles and Dipoles 
and Rhombic. The fields are calculated using a single 
point source for each element with length less than half a 
wavelength. For elements longer than this the number of 

2 0 point sources increases proportionately with the number of 

half wavelengths in the elements. 

As with the aperture antennas before the wire 
antennas can be properly modeled it is necessary to obtain 
the 3dB bandwidth in the two orthogonal radiation 
25 patterns. Once this has been done an estimate may be made 
of the number of elements N in the antenna. 

With the number of elements equaling: 

TV = — where 1 is the length of the antenna and X is 

A, 

wavelength 

3 0 the length and space of each of the wire elements can be 

determined. 

The method then requires a calculation of the 
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mutual coupling between all elements in the array and the 
assembly of the N x N impedance matrix. 

To obtain the current I in each element the 
inverse impedance Z matrix is multiplied by the column 
5 voltage V vector (0 for parasitic elements, 1 volt for the 
driven element) • 



h 

| | 

In 

10 

This gives an example for determining the current 
for an n element array. 

A Huygen's wavelet point source is then assigned 
15 to the current locations of each element in the antenna. 
The magnitude and phase of the point source is directly 
proportional to the current calculated previously. 

Assuming a sin (9) dependence for each wire element 
(9 is measured from the direction of the elements) , the 
20 contribution of each point source is added to every point 
in space in the vicinity of the antenna. This includes 
both near and far field affects. Thus the effect of each 
point source on a point in space is calculated in 
determining the overall power density level at the point 
25 in space. This is similar to the approach used in the 
Huygen's wavelet method used for the aperture antenna 
algorithm. 

Field strengths are also scaled dependent upon the 
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power supply to the antenna. 

The effect of metallic and dieletric support 
structures can be modeled by point sources using image 
theory. These images contribute to the array of point 

5 sources in the model and the field intensity calculations 
are made as before. 

In its preferred form the above described 
embodiments of the invention are implemented by computer 
software. The software is preferably configured to store 

0 data relating to different types of antenna. This data 
would include manufacturers antenna pattern files in a 
number of different standard formats. In this way it is 
possible to compare an antenna being modeled with existing 
data. 

5 If an antenna being modeled does not fit an 

existing manufacturers antenna pattern file the software 
is able to receive measured data relating to antenna 
patterns and create an appropriate file. This can be 
added to a main database for storing antenna pattern 

0 files. 

Typical antenna pattern properties which are 
stored include pattern type, frequency for that pattern, 
system loss, resolution, linear averaged, pattern cut, 
pattern type, electric tilt and effective gain. 

5 In addition to the above files storing basic 

antenna patterns may include data relating to the 
horizontal beam width, the vertical beam width and the 
front to back ratio. It is preferred that once the main 
characteristics of an antenna have been characterised the 

0 antenna can be modeled according to one of the above 

described methods. A power density graphical or visual 
representation then can be obtained and the software is 
able to highlight the near field and far field radiation 
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patterns using different colours. These differences are 
highlighted in black and white in Figures 3 to 6 • 

It is preferred that the software is able to 
select a number of resolution options so that the 
5 resolution of the power density patterns can be varied. 

It is also preferred that the software is able t 
record safe power density levels and restrict the images 
produced for near field and far field radiation to the 
safety values- This diagrams can be produced which show 
10 only the hazardous levels of radiation. 
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CLAIMS 

1. A method for determining £ield radiation 
levels for a radiating device comprising determining far 
field radiatiofi characteristics of a radiating device, 
5 providing a mofLel of the radiating device, which model 
approximates tfie determined far field radiation 
characteristic^ and determining a near field radiation, 
characteristic from the model for at least one point in 
space.. 

10 2 m The method as claimed in dlaim 1 including 

the step of determining a boundary between the near field 
and far field Radiation of the radiating device. 

3. The method as claimed in claim 1 or 2 
including the £tep of determining near field radiation 

15 density from t|ie model. 

4. The method as claimed in claim 3 including 
the step of determining power density level over a 
plurality of positions in space. 

5. The method as claimed in claim 4 including 
20 the step of determining beam width characteristics of the 

radiating device in two orthogonal far field radiation 
patterns • 

6. The zttethod as claimed in claim 5 including the 

step of determining the 3dB beam width in two orthogonal 

25 far field radiation patterns.. 

7. The method as claimed in claim € including 
the step of determining physical characteristics of the 
radiating device to determine the far field radiation 
characteristic? . 

30 8. The method as claimed in claim 7 including 

the step of providing a model including Representing the 
device by a plurality of radiation sources. 

9. The method as claimed in any one of the 
preceding claims wherein the radiating device comprises a 

35 wire antenna. 

10. The method as claimed in glaim 9 including 
the step of providing a model including Representing the 
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radiating device by a plurality of wire elements. 

11. Tlie method as claimed in claim 10 including 
the step of estimating the length and spacing of each wire 
element forming the radiating device. 
5 12. The method as claimed in claim 11 wherein 

each wire element is represented as a radiation source. 

13. The method as claimed in claim 12 including 
the step of calculating mutual coupling between all the 
wire elements . 

10 14. The method as claimed in claim 13 including 

the step of assembling an N by N impedance matrix and 
calculating the voltage for each element to determine the 
current in each element . 

15. The method as claimed in claim 14 including 
15 the step of multiplying the inverse impedance matrix by 

the column voltage vector to determine the current in each 
element • 

16. The method as claimed in claim 15 including 
the step of assigning a Huygen's wavelet point source to 

2 0 each element and calculating the magnitude and phase of 

each wavelet point source from the current determined. 

17. The method as claimed in claim 16 including 
the step of summing the contribution of each point source 
to each point in space within the near field. 

25 18. The method as claimed in claim 1, 2 or 3 

including the step of providing a single point source for 
each element with a length less than half a wavelength. 

19 . The method as claimed in any one of claims 
1 to 7 or claim 17 wherein the radiating device is an 

3 0 aperture antenna. 

20. The device as claimed in claim 19 including 
the step of determining the physical characteristics of 
the radiating device and providing a model including 
representing the aperture by at least one Huygen's wavelet 

35 source. 

21. The method as claimed in claim 20 including 
the step of summing the contribution from each wavelet 
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source to each point in space, 

22. The method as claimed in claim 21 wherein 
the power density level at any point in space is 
determined using the formula 

&/+2.I5 

^ , PoweratAntenna * 1 0 10 



23. A method of estimating radiation power 
density of electromagnetic radiation comprising the steps 
of identifying a radiating device, representing the 

10 radiation device as a plurality of point sources which 
radiate electromagnetic radiation, estimating power 
density level at a plurality of positions in space for 
each point source and determining the total power density 
level at each position, by summing the contribution of 

15 each point source to the respective positions in space. 

24. The method as claimed in claim 23 including 
the step of displaying the power density level for a 
plurality of positions. 

25 . The method as claimed in claim 24 including 
2 0 summing the power density level determined at each 

position for all point sources representing the radiating 
device . 

26. The method as claimed in claim 25 including 
the step of calculating far field and near field tapering 

25 characteristics for each position. 

27. The method as claimed in claim 26 including 
the step of calculating the power density level at a point 
in space by using the power density formula 
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for far field radiation and modifying tlie far field power 
density formula for near field radiation, which 
modification affects the antenna gain, power sent to the 
antenna and the distance from the antenna to the point 
source . 

28. A method of determining field radiation 
levels for a radiating device comprising the steps of 
determining far field radiation characteristics of a 
radiating device, determining the boundary between near 
field and far field radiation, determining the 
displacement of a point in space relative to the closest 
point on the radiating device and calculating the power 
density level at the point in space. 

29. The method as claimed in claim 28 including 
the step of modeling the radiating device as a plurality 
of point sources . 

30. The method as claimed in claim 29 including 
the step of applying a closest point algorithm to 
determine the power density level at each point in space. 

31. The method as claimed in claim 30 wherein 
the closest point algorithm determines the displacement of 
the point in space from the closest point on the radiating 
device . 

5 32. The method as claimed in claim 31 wherein 

the closest point algorithm calculates X, Y, Z 
displacement vectors from the point in space to the 
closest point on the radiating device and calculates 
azimuth and elevation angles to the closest point. 

0 33. The method as claimed in claim 33 wherein . 

the closest point algorithm determines the orientation of 
the radiating device and scales the power density level 
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determined according to the orientation of the radiation 
device . 

34. The method as claimed in claim 34 wherein 
the closest point algorithm calculates the power density 
5 level using the power density formula and incorporates any 
modification factor applicable if the point in space is in 
the near field. 
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